Magnetic and Pairing Correlations in Correlated Electron 

Systems 



Gang Su* 

Institut fiir Theoretische Physik, Universitdt zu Koln 
Zillpicher Strasse 77, D-50937 Koln, Germany 

Abstract 



Magnetic and superconducting pairing correlation functions in a general 
class of Hubbard models, the t-J model and a single-band Hubbard model with 
additional bond-charge interaction are investigated, respectively. Some rigor- 
ous upper bounds of the corresponding correlation functions are obtained. It 
is found that the decay of the spin-spin correlation function with temperature 
in the general Hubbard models can not be slower than the squared inverse 
law at low temperatures and the inverse law at high temperatures, while the 
on-site pairing correlation function can not be slower than the inverse law. 
An upper bound for the average energy of the t-J model is found. The up- 
per bounds for the spin-spin and the electron pairing correlation functions in 
the t-J model as well as in the Hubbard model with bond-charge interaction 
are also obtained. These bounds are expected to provide certain standards 
for approximate methods. In some special cases, these bounds rule out the 

possibility of corresponding magnetic and pairing long-range order. 
PACS numbers: 71.27.+a, 75.10.Lp, 74.20.-z 



1 



Typeset using REVT^ 



I. INTRODUCTION 



The study of superconductivity and magnetism in strongly correlated electron systems 
has been receiving intense interest recently. This may be attributed to the discovery of 
high temperature superconductivity Since the idea to explain superconductivity in the 
framework of strongly correlated electron systems (Hubbard model and its variants) was 
proposed [1], numerous investigations on these systems have been done. For the many-body 
problems are difficult to solve, only a few rigorous results are known in literature [2]. Most 
of people used approximate or numerical methods to investigate spin and pairing correlation 
functions, and usually achieved some different, sometimes even conflict results. In addition, 
some difficulties in numerical methods remain unsolved. Therefore, a variety of interesting 
questions in these systems remain ambiguous to date. In this situation, it is really necessary 
to search for some exact results which, on the one hand, can be used to examine the validity 
of some kinds of approximations or numerical calculations, and on the other hand, can help 
us to further understand physical properties of these systems. 

In this paper, we shall study magnetic and pairing correlation functions in the following 
systems: a general class of Hubbard models, the t-J model and a single-band Hubbard 
model with additional bond-charge interaction. By using Bogoliubov's inequality, we give 
some upper bounds for the studied correlation functions, which, at some extent, may provide 
certain checks and standards for approximate methods. In some special cases, these bounds 
rule out the possibility of the corresponding magnetic long-range order (LRO) and the on-site 
pairing LRO. 

The rest of the paper is organized as follows. In Sec. II a general class of Hubbard models 
are studied. Symmetries of correlation functions are discussed. The upper bounds for the 
magnetic correlation function and the on-site pairing correlation function are presented, 
respectively. In Sec. HI and IV, the magnetic and pairing correlation functions in the t- 
J model and a single-band Hubbard model with bond-charge interaction, are discussed, 
respectively. A simple upper bound for the average energy of the t-J model is obtained. The 
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upper bounds for the spin-spin correlation functions and the pairing correlation functions in 
the two models are given. In Sec.V a summary of the results are presented. 

II. HUBBARD MODELS 

Consider a general class of Hubbard models on a d-dimensional lattice with M (even) 
sites. The Hamiltonian is 

H = ^ tj.,r' (^r^r-' + ^JV) + X! Uralttrblbr — ^/ir(ajar + &J&r), (2.1) 
r,r' r r 

where try is the hopping matrix element, and satisfies t*^, — tr\r- «r (br) annihilates a 
spin up (down) electron at site r, Ur is the local spin-independent Coulomb potential, and 
/ir is the position-dependent chemical potential. No other a priori assumption (apart from 
those indicated explicitly in the context) is needed. This is thus a very general form of the 
Hubbard model. 

Define the local spin operators as follows: 5"+ = ajft^, S~ — feja^, — |(n" — n^) 
with < = alar and = blbr, and [S+,Sr>] = 2S^6ry, [S^^S^,] = The global 

spin operators 5^ = X)r ^^'^ = Y,r S,- obey the usual SU(2) symmetry and satisfy 
[S^, H] = [S', H] = [S^, H] = 0, where = 1{S+S- + S-S+) + {S'f has the eigenvalue 
S{S + 1). The particle number N — X^^(ajar -l- blbr) — + is conserved, and commutes 
with the Hamiltonian. Define the rj operators as follows: 77+ = ajfej., rj~ = brUr, rj^ = 
^{rir — 1) with rir — + n^, and [r]:^,r]~,] — 2r]r5ry, [^^,^r'] — ^'ut^ry- Below we will 
investigate the spin and the on-site pairing correlation functions of the model, and these 
definitions are necessary for subsequent analyses. 

A. Symmetries of correlation functions 

First let us for convenience write down three well-known unitary operators explicitly [3] , 
which were frequently cited in hterature, and were usually applied to study the trans- 
formed systems connected by them, but their explicit forms are not obvious. The op- 
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erator Uo — Yiri^r — e(r)6j) with e(r) = (—1)'' and UqUo — 1 designs the well-known 
particle-hole transformation [4]: UoaMo = cLr, UohrUl = e(r)6j, which makes Ur —Ur 
in the Hamiltonian (2.1) with /i^ = ^ if = ~t for nearest neighbors and zero 
for others, hke the standard single-band Hubbard model. Here (—1)'' can be under- 
stood as a factor e**"^'' with Q = {n,n,...) in two or three dimensions. The operator 
^1 = rirlctr — cil){br — with ulUi — 1 designs another symmetric particle- hole trans- 
formation: UiajLll — a\, UihrUl — h\ which makes — > —S~, — > —S^, 77^ —rj~ 
and tr,r' —tr',r in (2.1) with additional constants and proper adjustment of the chemical 
potential. The operator U2 = exp[YJ2ri^l — bl)i^r — br)] with U2U2 = 1 exchanges the 
spins (up-down symmetry): U2a\U2 — b\, which leaves the Hamiltonian (2.1) unchanged, 
but makes S~ and —S^. Of course, these unitary operators can be used either 

individually or in a combined way. 

We now give some relations for thermal correlation functions. It is noteworthy that 
some of which are trivial, but some are less obvious and thus worth to write down. As the 
expectation value of the commutator, [n", H] = J2r'{'tr,r'<4-'^r' ~ tr'.rCt^Qr)) vanishes, we get 

'^try < alttr' >= ^ V,r < Cll'Ctr >, (2-2) 

where < • • • > denotes the thermal average. Using the up-down symmetry, connected by 
U2, we have 

"^try < blbr' >— ^ V,r- < ol'dr > ■ (2.3) 

Note that (2.2) and (2.3) are general, not limited to the translation invariant case, for we 
have not made use of any spatial symmetry of the lattice. The expectation value of the 
commutator [S+S^, S'] = 2S^,S^ - S+S^ yields 

< S+,S- >= 2 < S',,S', > (2.4) 

due to commuting with H. This symmetry is essential, for it is local, and is valid for 
any r and r'. Similarly, we can obtain a lot of such symmetries, for instance, 

< S^^S^ >= 0, < >=< >= 0, < ^X' >= 0, 
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< >= 2 < >, 

< n°,n^ >=< n"n(;, >, < a,.6r >=< 6^4 >= 0) 

< (-5,^)'"*+' >= (2.5) 

with m — 1,2, 3,..., etc. It can be seen that some of above relations are obvious for transla- 
tion invariant case, but less obvious for the system without the translation invariance, which 
are probably useful for numerical calculations, and meanwhile give some restrictions on ap- 
proximate methods. Moreover, one may apply the unitary transformations, as mentioned 
above, to Eqs.(2.5), and can obtain the corresponding symmetries of correlation functions 
for the transformed systems. We will apply them below. 

B. Magnetic correlation 

Let us study the transverse spin correlation function < S^,S~ > for r ^ r', which is 
related to the longitudinal correlation function through (2.4), by means of Bogohubov's 
inequality [5] 

I < [^, S] > 1^ < ^ < {A, At} X [[B, H], St] >, (2.6) 

with P — {ks — 1) the inverse temperature. Note that the relation between spin 
correlation function and magnetic LRO has been extensively established thirty years ago [6]. 
Since 

we have 

< HI S',] >= - E V < (^W >> 0, (2.7) 

where the non-negativity of (2.7) comes from the fact that the inner product {B, B*) > 0[5,7]. 
By the Schwartz inequahty | < ^It^ > | < V< A^A >< B^B >, we observe that 

<{S+S;,S+S;,}><8 (r'j^r). (2.8) 
5 



To obtain inequality (2.8) one has to substitute the definitions of spin operators into the left- 
hand side of (2.8), and then apply the Schwartz inequality to electron operators repeatedly 
by noting that | < a|!,ar > | < 1 and \ < bl,br > \ < I. Therefore, setting A — S^,S~ [r ^ r') 
and B — in (2.6) and noticing (2.7) and (2.8), we get a bound 

I < S^S; >f< -4/3 try < alar' > (2.9) 

r'i^r) 

for r' r. One may note that the index r' in the right-hand side (RHS) of (2.9) has been 
eliminated due to using the Schwartz inequality in (2.8). The same situation occurs in the 
following. Obviously, as tr',r{r' 7^ r) — > 0, then | < 3^,3^ > | = 0, which implies that 
there is no spin-spin correlation in the atomic limit, i.e., no magnetic (ferromagnetic and 
antiferromagnetic) LRO occurs in this case. One may notice that the RHS of inequahty (2.9) 
depends only on off-diagonal correlation function < aja^' > and hopping matrix element try, 
independent of local Coulomb potential Ur- 
Define operators 

t - - - t - t 

air' ^ ^r',r-^r' " ^^.r'^r, Oiry = t^yOl, - t^, ^^al- (2-10) 

Then < al^ary >> gives 

try < alttr' >> - ^ |V,r|(< alittr' > + < aja^ > —2). (2-11) 

Substituting (2.11) into (2.9) we have 

I < S:^,Sr > 1^ < 2/3 l*r',r|(2- < alar' > - < alar >). (2.12) 

If the system has translation invariance, then we from (2.12) verify rigorously a trivial fact 
I < 8:^,3- > I = for r' 7^ r at full-filling. Although the bound (2.9) is lower than (2.12), 
sometimes the latter is also expected to be useful. 

Since the RHS of (2.9) is intimately related to < alar' >, we now discuss it. For r' ^r, 
we have [ajor', aja^] = — aja^'. By setting A = alar' and B — alar in (2.6) one obtains 

I < alar' > 1^ < -/3 < « - </)^ > Y try < alar' >< -2/3 ^ try < alar' > 

r'{^r) r'{^r) 
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for r' 7^ r. This is a recursion inequahty. It turns out to be 

I < alar' >\<2(3 I Vl- (2-13) 
Moreover, we are also able to obtain a bound 

I < alar' > I < Y^(< < > - < «, >)(!- < n« >) (2.14) 

for r 7^ r'. The two bounds combining (2.9) may help us to understand more about the spin 
correlation function. 

Let us turn to a special case for the moment. Assume that = —t with t > for 
r, r' being nearest neighbors, and otherwise, like the standard single-band Hubbard model 
but including the local Coulomb potential Ur- We introduce the Fourier transform of a J as 
al = Z^p OpC"*^'', where the summation on p runs over the dual lattice defined by the 
boundary conditions. By summing over r r') on both sides of (2.9) and inserting the 
Fourier transform into it we obtain Z)r(^r') I < S^iS~ > P < 2t/3J2p,s < Up > e'^"^, where 
\d\ denotes the lattice spacing between nearest neighbors, and = a^a^ + 6^6^. One may 
observe that when < >= 1 or constant, we have | < > | = for r 7^ r' due to 

J2p e*^'' — 0. In other words, the single-band Hubbard model with local Coulomb potential 
does not exhibit magnetic LRO at finite temperatures at < >= 1 or constant. This result 
is independent of the sign of Ur and for arbitrary dimensions. Although the condition of 
< Up >— 1 or constant is very special, we rigorously rule out the possibility of magnetic 
LRO in the case. 

If we donot bound < {S^,S~,S:^S~,} > by the Schwartz inequality in (2.8), we have an 
expression 

< ><-< 

where we have used (2.5) and < {A^, A] >> 0. Then from (2.6) we get 

I < 5+5- > I < /3Q{r) + ^jQ{ryp^ + pP{r,r')Q{r) 
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with Q{r) = - Er'(^r) try < alar' > and P{r, r') = | < S^^S^^S+S' > | < 1 for r r'. By 
noticing (2.13) we can obtain an upper bound for the spin-spin correlation function 

I < S^,S; > I < G{f3) = 2/3^ R{r) + ^ Aj3^R{rY + 2j3^R{r), (2.16) 

with i?(r) = (I]r'(7^r) l^r,r'|)^- To show the bound unambiguously, we assume that — i 
for r, r' being nearest neighbors and otherwise. We plot G(/3) versus temperature as 
shown in Fig.l, where (5"^ is in units of 2i, and the coordinate numbers are taken as 6, 4, 2, 
respectively. Since | < S'^S'" > | < 1, we only plot the interesting part. Prom Fig.l one 
may see that the bound decreases rapidly with increasing temperature. When /?~^ > 20, the 
bound decreases slowly, eventually to zero, as temperature increases. Evidently, the decay 
of the spin-spin correlation function with temperature can not be slower than the squared 
inverse law at low temperatures and the inverse law at high temperatures. Although the 
bound as well as the bound (2.12) can not give a general reply if the system can exhibit 
magnetic LRO or not, it may shed useful light on examining the validity of some kinds of 
approximations and numerical results, especially on the dependence of temperature of spin- 
spin correlation functions. We note that the spin-spin correlation function of the single-band 
Hubbard model was studied numerically on small sizes (4 x 4) of a square lattice at half- 
filling at low temperatures[8]. Their calculated results are found to be smaller than the 
present bound, as indicated in Fig.l. The reason for this discrepancy may be that in spite of 
the finite-size effects in numerical calculation, the present bound is suitable for macroscopic 
sizes of lattices and is better for high temperatures. 



C. Pairing correlation 

To investigate the on-site superconducting correlation, we need to calculate the on-site 
pairing correlator < ?7^?7^ >=< aj6j6r'flr' > in the off-diagonal long-range (ODLR) limit[9] 
|r — r'l — > oo, namely, off-diagonal long-range order (ODLRO)[10]. As before, we use the 
Bogohubov's inequahty. Choosing A — rj^^rj', and B — rj^ in (2.6), and noticing that 
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< H],r)::] - E tr',r < alar >, (2.17) 

we have the bound 

I < VrVr' >\^ <-/3 ^r',r < a^Qr >, (2.18) 

where we have used the Schwartz inequahty to bound < {'r}rVr'^Vr'''lr} ^— 2- One may 
observe that if try — > 0, then | < ry+r^^ > | — > 0. This suggests that no on-site pairing 
correlation in the general Hubbard model exists in the atomic limit. 
Since 

E < al,ar >= X] < alttr > -N^ = M < > -N^, (2.19) 

r,r'{r^r') r,r' 

with < Uq >—< alao > the number density with zero momentum of spin-up electrons, and 
further assuming that tr'^r = t — const., from (2.18) one obtains 

^ E I < > r < - y(< no > -p), (2.20) 

with < no >—< alao + blbo > and p — J2r < ^r > /M. We note that if t > and < n^ >> p 
or t < and < no >< p, then | < aj6j6r'ar' > | = for r 7^ r' from (2.20). In other 
words, the system can not exhibit the on-site pairing condensation in the aforementioned 
circumstances. Prom the derivation, one may note that the electron hopping terms plays a 
key role in pairing condensation phenomena in itinerant electron systems. Besides, one may 
observe that the sign of t also has the effect on the final result, as shown above. Of course, 
this argument can also apply to (2.9). 

If we exchange A and B in the derivation of (2.18), then we can get | < 77^77^ > P < 
13 < {r^^f >< [[B, H],B^ >. One may see that if 

< Knl >< ^(< Ur > -1), (2.21) 

it gives I < 77+77^ > I = 0. Namely, under the condition of (2.21), the Hubbard model can 
not appear the on-site pairing condensation at finite temperatures. The condition (2.21) is 
not pecuhar, e.g., the case of half- filling with < n^n^ >— falls in it. 

9 



If tr,r' — —t for r,r' being nearest neighbors and/or next nearest neighbors and for 
others, then similar to discussions for spin correlation functions we also have | < aj6j6j./aj./ > 
I = for r 7^ r' at < rip >= 1. We notice that Veilleux et al[ll] have recently studied the pair 
correlations of the Hubbard model with next nearest neighbor hopping by using quantum 
Monto Carlo method. Their consequences are qualitatively in agreement with the present 
rigorous result in their studied parameter region. 

By substituting (2.13) into (2.18), we get a rigorous upper bound for the off-diagonal 
element of the two-particle reduced density matrix 

I < alblbr'ttr' > I < V2P \*r',r\: (2.22) 

which should also be valid in the hmit \r — r'\ ^ oo. This bound suggests that the decay 
of ODLRO with temperature is not slower than the inverse law. Although the bound is too 
loose to compare with the numerical data on small clusters at low temperatures [8], it gives a 
hint for the pairing correlation function in the thermodynamic limit. Since < aj^jfer'^r' > is 
related to the superconducting order parameter < Urbr > (quasi-average) by the asymptotic 
property[12] < alblbr'ar' >— > | < ttrbr > P in the ODLR limit, (2.22) provides a standard 
for approximants in calculating the temperature-dependence of the pairing order parameter. 
Here we would like to point out that one may obtain the similar bounds as (2.22) for other 
pairings, for instance, the extended s-wave, d-wave, and so forth[13]. 

III. T-J MODEL 

This model has been extensively studied in recent years, but the rigorous result is rare, 
except that the one-dimensional (ID) supersymmetric model (J = ±2t) is exactly solved 
using Bethe ansatz[14]. Many approximate or numerical results on magnetic and pairing 
correlations in high dimensions in this model are quite different so far [8]. We study the 
following Hamiltonian 

Ht-J = X] (4«r-' + blbr') + J X] i^r- Sr' - ^nrUr') (3.1) 

<r,r'> <r,r'> 
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on a d-dimensional lattice, where the notions are the same as in Sec.II, < r,r' > are nearest 
neighbors, J > (we here consider J, without loss of the generality, as an independent 
parameter) and t > 0. The different forms of the model has been discussed elsewhere[15]. In 
this model, we assume that the double occupancy of every site is excluded. In other words, 
each lattice site is constrained to have either one electron (with spin up or down) or none, 
as usual. It can be seen that the system has SU(2) spin symmetry. We may also obtain 
some symmetries of correlation functions as in Sec.II. In this section we will first derive 
an upper bound for the average energy, then study the spin-spin correlation function, and 
finally discuss the nearest neighbor pairing correlation function. 



A. Upper bound for the average energy 

Prom (3.1) we find 

< [[S+, Ht-j],S;] >= t ^ < al,ar + bX' > -2J ^ < 2S^S^, + S+S; >, (3.2) 

r' r' 

' <r> ' <r> 

where r'^r> denotes the summation on r' running over nearest neighbors of r. (3.2) then 
implies 

J2 < S^'S; >< E < + blbr' >, (3.3) 

where we have used (2.4) and the non-negative property [5, 7] of (3.2). By noticing (3.3) and 

< S+S~, >=< S:^S~ > one gets < Ht-j >< -ft Y.s,p < rip > e'^^ -| J2<r,r'> < nrUr' > . At 
temperature T, we have the average energy (internal energy) Eq =< Ht-j >. On the other 
hand, the non-negativity of < [[a^fcr, -f^t-j], ^'JaJ] > yields —tY,<r,r'> < a\.ar' -\- b\.br' >< 
■f I]<r,r'> < nrUr' > —^Nz with z the coordinate number, and thus Eq < jgJ2<r,r'> < 
rirUr' > —^Nz. Furthermore, since < (n^ — rir'Y >> and noting that < n"n^ >= due 
to the restriction of no doubly occupied sites, one has < n^rir' >< |(< rir > + < rir' >)■ 
Substituting these facts into Eq we have 

Eo < -\jzN. (3.4) 
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We like to point out that the bound (3.4) is generic, not hmited to the translation invariant 
system, and is valid for arbitrary filling fraction and arbitrary dimensions. If the system 
has the singlet ground state, then Eq ai T ^ can be regarded as the ground state energy, 
and < ... > thus means the average in the ground state. If the ground state of the system 
is degenerate, Eq can also be understood as the ground-state energy for all ground states. 
At half-filling, the t-J model reduces to the Heisenberg antiferromagnetic model. In ID 
the ground state energy is well known to be Eq/M — — JZn2[14], which clearly satisfies 
(3.4): Eo/M = -0.693147 J < -0.5 J. Away from the half-fiUing, the ground state energy 
of the supersymmetric t-J model also comply the bound (3.4), as shown in Bares et al's 
paper[14]. In 2D, the numerical result for estimating the ground state energy of one hole 
in the interval 0.2 < J/t < 1.0 on small clusters (4 x 4) is Eq/M = -3.17t + 2.83t{J/tf-^^ 
[16]. However, their Hamiltonian does not contain the term J2<r,r'> nrUr'- If this term is 
taken into account, their numerical results would be in agreement with the present bound. 
Since this bound is a rigorous result, it can be expected to provide checks for approximate 
and numerical methods, especially in high dimensions. 



B. Magnetic correlation 

The magnetic structure factor m{p) is given by 

^ip) -j^H ^''^''"''^ < s^si, >=< s;s% >, (3.5) 

r,r' 

where 5^ = ^ -^^e*?"". By symmetry, < ^^^Ip >= \ < S+S- >. Prom (3.3) we obtain 

E[^(P) - ^ < >]7f ^ (3-6) 

with 7p = I^^e*^''. Note that m{p) > 0. Inequality (3.6) gives a severe restriction on m(p) 
in the t-J model. If < Up >= 1 or t = 0, we have J2p'>TT'ip)'lp < 0. Since the existence of Neel 
order corresponds to m(p) containing a 5 function at Q = (tt, tt, ...) in the infinite-volume 
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limit[17], and let be the coefficient of tfiis delta function, we from (3.6) get a bound 



t 



m 



(3.7) 



P P^Q 

It has been shown that the 3D half-filled t-J model, i.e., Heisenberg antiferromagnetic model, 
has LR0[17]. Away from half-filling, (3.7) may shed some light on the antiferromagnetic 
order of the t-J model. For a square lattice, 7(5 = —4. Then we have 



t 



m 



P^Q — ' P 

If we obtain a bound for Y.p^q ^{p)lp-i then we can say something about the antiferromag- 
netic LRO in the model, which will be left for future study. 

We choose A = S-,S^ and B = S+ in (2.6). Then [A,B\ = S;,S^^ for r 7^ r' . Since 
I < {S~iS^,S^S^i} >< I by the Schwartz inequality, from (2.6) and (3.2) we have 



< s+s; > I < 



f (t V < al,ar + hlhr' > -4J V < S+S- >) 



<A^[t + 2J^^{t + 2J^'-^{t + 2 J^))). 



(3.8) 



This inequality gives an upper bound for the spin-spin correlation function in the t-J model. 
Particularly, as t — 0, the model reduces to the Heisenberg antiferromagnetic model, and 
(3.8) becomes 

I < S^S; >\<^z. (3.9) 

That is, the temperature-dependence of the spin-spin correlation function can not be slower 
than the inverse law in the Heisenberg antiferromagnetic model. 



C. Pairing correlation 

Now we come to discuss the pairing correlation function. Since there is no doubly 
occupied sites in the system, the on-site pairing correlation should be vanishing. We in 
the following cosider the nearest neighbor pairing order parameter < a^^ri > with ri,r2 
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being nearest neighbor sites. For this purpose, we have to add a U(l) symmetry breaking 
term H' — —a J2<ri,r2>i'^ribr2 + ^Ja'^n) ^ ^^^^ ^ chemical potential term — /x J2r into 
the Hamiltonian Ht-j. Let B = ar^br2, and A = aj^a^i + b^^br^- Then [A, -B] = —2arJ)r2i < 
[[ur^br^, Ht-j],bl^alJ >< D{ri,r2) with D{ri,r2) = ^(Er<.i> < 4i«r- > +Er<.2> < bl^br > 
)-(J + 2/x)(l- > - << >) + '2J K^y^^y +6Jz-a{Er,^^,albl+Er,r^,blalJ, 

and < {^4,74''^} >< 8. Substituting these results into (2.6) we obtain the bound 

\<ar,br2>\^ <PD{n,r2) (3.10) 

as q; — > in the thermodynamic limit, where we have used the Schwartz inequality to 
bound those terms with four creation and annihilation operators. One may note that < 
0'ribr2 > in (3.10) is the Bogoliubov's quasi-average[12]. On the other hand, according 
to Bogoliubov's argument [12], the off-diagonal element of the two-particle reduced density 
matrix < 6|!^a|!, 0^16^2 > has the asymptotic behavior 

< bl,^al,^aribr2 >^< >< a^br^ > (3-11) 

in the ODLR limit[13] \{r2,r[) — (r2,ri)| — > 00. It is worth mentioning that (3.11) is not 
incompatible with Haag's spatial cluster theorem[19]. In the translation invariant system, 
the temperature-dependence of ODLRO should thus obey 

><PD{n,r2), (3.12) 

i.e., the decay of the ODLRO for nearest-neighbor pairs with temperature in the t-J model 
with translation invariance can not be slower than the inverse law. This bound thus offers 
a check for some approximate results on the temperature-dependence of superconducting 
order parameter. 

IV. HUBBARD MODEL WITH BOND-CHARGE INTERACTION 

The Hamiltonian of the model is given by 

i?6_c = X] i^l^r' + blbr' + h.C.) + ajorfcj&r- 

<r,r'> r 
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+X [{alar> + h.c.){nl + n^,) + {hlW' + /i.c.)« + <0] - A* + ^J^), (4.1) 

<r,r'> 

where X is the bond-charge interaction, and other notations are standard, as usual. This 
model has been extensively discussed by Hirsch[20], Bariev et al[21], de Boer et al [22] 
and Schadschneider[23]. It has been shown[22] that the rj pairing state with ODLRO is the 
ground state if ^7 < — 2Z|t|, and can be solved exactly in one dimension[21-23], &st — X and 
jjL — U/2. For certain values of X and large densities of electrons (small doping) the bond- 
charge interaction may lead to an attractive effective interaction between the holes within 
the framework of BCS mean- field theory[20]. In this section, we will give a few rigorous 
bounds for the spin-spin and the on-site pairing correlation functions of this system. 

Like the way exploited in preceding sections, it is not difficult to obtain the following 
two bounds 

I < S^S;, > I < [8(3{t J2 < 4'(^r >-Xj2< alMnr + n^') >)]^ (4.2) 

r' r' 

' <r> ' <r> 

and 

I < alblbr'ttr' > I < [2P{t J2 < > -X J2 < al,ar{nl + nl,) >)]^ (4.3) 

r' r' 

' <r> ' <r> 

for r 7^ r' by setting A — S^S~, and 77^77^ with r ^ r' and B — and 77^ in (2.6), 

respectively, where we have used the spin- flip symmetry (W2). Evidently, if X satisfies the 
following condition 

X Yl < ol'dririr + n^O >>tY< al,ar >, (4.4) 

r' r' 

' <r> ' <r> 

then there is no spin-spin correlation and on-site pairing correlation. To assure the existence 
of magnetism and superconductivity in the system i^b-o the condition in opposite direction 
of inequality (4.4) must be hold, which gives a restriction on the values of X. 

For a special case t = X, the system possesses the symmetric particle- hole symmetry, 
connected by the unitary operator Ui, as discussed in Ref.[23], the 77 pairing symmetry [22] 
at half-filling, and so forth. This model is really interesting, and the details will be presented 
in a seperate publication. 
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V. SUMMARY 



We have rigorously investigated magnetic and superconducting pairing correlation func- 
tions in a general class of Hubbard models, the t-J model and a single-band Hubbard model 
with additional bond-charge interaction by means of the Bogoliubov's inequality, respec- 
tively. Some corresponding upper bounds are obtained, which are expected to provide 
certain checks and standards for approximate methods. In some special cases, these bounds 
rule out the possibility of corresponding magnetic and pairing LRO. 

For the Hubbard models, we obtained an upper bound for the spin-spin correlation 
function, which indicates that the decay of the correlation function with temperature can 
not be slower than the squared inverse law at low temperatures and the inverse law at 
high temperatures. Prom these bounds we observe that there is no magnetic LRO in the 
atomic limit and in the case with the momentum distribution function being constants for 
the model with only nearest neighbor hoppings. An upper bound for the on-site pairing 
correlation function was obtained, which suggests that the decay of the on-site pairing 
ODLRO with temperature is not slower than the inverse law. Since our method is rigorous, 
the present result may be applied to clarify some contradication in approximate calculations. 
In addition, we found that there is no on-site pairing correlation in the atomic limit and in 
the cases either t > and < Uq >> p or t < and < no >< p (see (2.20)) or < > being 
constants for the single-band model but with local Coulomb interaction. It is emphasized 
that all obtained bounds are independent of the local on-site Coulomb interaction and are 
valid for arbitrary dimensions. 

For the t-J model, we obtained an upper bound for the average energy (internal energy) 
for arbitrary electron fillings. Whenever the ground state of the system is unique or not, the 
upper bound at T ^ can be regarded as that of the ground-state energy. Since the bound 
is rigorous, it provides a standard for approximate and numerical methods. We also obtained 
a lower bound for the Neel order, which may shed useful light on the antiferromagnetic order 
of the system. An upper bound for the spin-spin correlation function was derived, which 
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implies that the decay of it with temperature in the model is not slower than the law away 
from half-filling and the inverse law at half-filling. An upper bound for the nearest-neighbor 
pairing correlation was obtained for the translation invariant system, which suggests that 
the decay of ODLRO with temperature can not be slower than the inverse law. The results 
hold for arbitrary dimensions. 

For the Hubbard model with bond-charge interaction, we obtained two bounds for the 
spin-spin correlation function and the on-site pairing correlation function, which gives a 
severe restriction on the values of the bond-charge interaction. 
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Figure Caption 

Fig.l The temperature dependence of the bound (2.16), where /?~^ is in units of 2t, and 
coordinate numbers are taken as 6, 4, 2 respectively, as indicated. The numerical data are 
taken from Ref.[8]. 
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